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h. 

1 
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mass  transfer  parameter  m^/p^  Ug 
constants  in  asymptotic  solution,  Eq.  7 

basic  flow  skin  friction  and  its  disturbance,  respectively 
p.'  /po.(i  = 1,2,3) 

Fourier  transform  of  h^ 

Wave  number  in  Fourier  transform 
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AP 

Q 

Re6,L 


T 

u,  v 

x,y 

B 

y 

6, 61 

r\ 


p,  v 


v 

P 

e1 

T 


mass  flow  rate  across  wall  surface  (positive  for  blowing) 
Mach  number 

undisturbed  and  disturbance  static  pressures,  respectively 
static  pressure  jump  across  incident  shock 
unit  solution  of  pressure  disturbance  Eq.  (3) 

Reynolds  number  based  on  conditions  ahead  of  shock  and 

either  boundary  layer  thickness  or  running  length  L, 

respectively 

absolute  temperature 

x and  y direction  velocity  components 

streamwise  and  normal  coordinate  distances,  respectively 
\l 1 or  V l-M^ 
specific  heat  ratio 

boundary  layer  and  displacement  thickness,  respectively 

y/e 

coefficients  of  dynamic  and  kinematic  viscosity,  respec- 
tively 

viscosity  temperature  dependence  exponent  ~ T**3) 
Fourier  transform  of  pressure  (Eq.  3) 
density 

boundary  layer  momentum  thickness 
shear  stress 


Subscripts 


Inviscid  properties  at  boundary  layer  edge 

i = 1,2,3  denotes  various  disturbance  regions  (Fig.  1) 

denotes  undisturbed  (not  stagnation)  flow  property 

Lighthill  friction  sublayer 

Eckert  reference  temperature  condition 

property  evaluated  at  wall  surface 
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TRANSONIC  SHOCK  - TURBULENT  BOUNDARY  LAYER 
INTERACTION  WITH  SUCTION  AND  BLOWING 

G.  R.  Inger  and  S.  Zee* 

Virginia  Polytechnic  Institute  and  State  University  , 

Blacksburg,  Va.  , USA 

Abstract 

? A basic  theory  of  weak  normal  shock  - turbulent  boundary  layer  inter- 
actions is  given  for  two-dimensional  non-separating  flows  including  mass 
transfer  across  the  wall  throughout  the  interaction  region.  Even  small 
amounts  of  suction^T-  mw/pe  Ue  ~ 1 to  ~5  x 10  ~^lare  found  to  signifi- 
cantly reduce  both  the  streamwise  scale  and  thickening  effect  of  the 
interaction  and  delay  the  onset  of  separation.  This  is  shown  to  be  a 
consequence  of  the  large  mass  transfer  effect  on  the  shape  of  the  in- 
coming boundary  layer  Mach  number  profile  away  from  the  wall.  Para- 
metric study  results  showing  the  influence  of  Reynolds  and  shock 
Mach  number  as  well  as  mass  transfer  parameter  on  the  interaction, 
plus  favorable  comparisons  with  various  experimental  data,  are  also 
presented. 


* 

Professor  of  Aerospace  Engineering  and  Graduate  Research  Assistant, 
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1.  Introduction 


The  study  of  transonic  shock-turbulent  boundary  layer  interactions  is  of 
well-established  importance  in  tne  aerodynamic  design  of  high  speed 
aircraft  wings,  transonic  wind  tunnels,  cascades  in  turbomachinery  and 
airbreathing  engine  inlets.  Consequently,  the  control  and  suppression  of 
interaction  effects  in  these  applications  by  suction  or  blowing  is  also  of 
great  interest  since  the  technology  of  boundary  layer  control  (BLC)  by 
surface  mass  transfer  has  advanced  to  a very  practical  status.  Although 
there  have  been  some  experimental  investigations  of  BLC  applied  to 
shock-boundary  layer  interaction  suppression  in  supersonic  inlets  * 
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transonic  wind  tunnel  side  walls  and  on  transonic  airfoils  ’ which 
have  established  the  desirability  of  the  idea,  little  systematic  basic 
study  has  been  done  in  the  transonic  regime , especially  with  regard  to 
establishing  a sound  theoretical  framework  for  the  problem.  It  is  the 
purpose  of  the  present  paper  to  develope  such  a fundamental  theory  for 
the  restricted  but  important  case  of  non-separating  flow  with  normal 
wall  (unvectored)  mass  transfer. 

Our  approach  is  based  on  extending  a recently-developed  analytical 
theory  of  weak  normal  shock  - turbulent  boundary  layer  interactions 
on  the  premise  that,  notwithstanding  the  existance  of  powerful  nume- 
rical methods,  there  will  be  a continuing  need  for  analytical  methods 
which  delineate  the  essential  physical  features  and  parametric  trends  of 
transonic  shock -boundary  layer  interactions.  In  Section  2 the  basic 
formulation  and  features  of  the  theoretical  model  including  surface  mass 
transfer  effects  are  given.  Section  3 presents  typical  numerical  results 
for  zero  mass  transfer  that  include  a number  of  heretofore  - unpub- 
lished Reynolds  number  effect  results  and  comparisons  with  experiment, 
followed  by  illustrations  of  the  mass  transfer  effect  on  important  inter- 
action properties  with  some  favorable  qualitative  comparisons  with  the 
little  data  available.  Section  4 concludes  with  a discussion  of  the  limi- 
tations of  the  theory  and  recommendations  for  further  studies. 
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2.  Theoretical  Formulation 


2.  1 Basic  Features  of  the  Interaction  Flow  Model 


It  is  well-known  experimentally  ’ that  when  separation  occurs,  the 
disturbance  flow  pattern  associated  with  normal  shock -boundary  layer 
interaction  is  a very  complicated  one  involving  a bifurcated  shock  pat- 
tern (see  Fig.  la),  whereas  the  unseparated  case  pertaining  to  turbulent 
boundary  layers  up  to  < 1.3  has  instead  a much  simpler  type  of 
interaction  pattern  which  is  more  amenable  to  analytical  treat- 

ment (Fig.  lb).  With  the  use  of  some  judicious  yet  physically- sound 
simplifications,  it  is  possible  to  construct  an  approximate  analytical 
theory  of  the  problem.  For  the  sake  of  orientation  and  completeness, 
a brief  summary  of  this  theory  will  now  be  given. 

The  flow  is  taken  to  consist  of  a known  unseparated  turbulent  boundary 
layer  profile  M (y)  subjected  to  small  transonic  disturbances  due  to  an 
impinging  weak  normal  shock.  Our  theoretical  model  of  this  interaction 
is  a simplification  of  the  small  disturbance  flow  structure  emerging 
from  an  asymptotic  analysis  of  the  compressible  Navier-Stokes  equations 

g 

at  high  Reynolds  numbers  , giving  a linearized  boundary  value  problem 
surrounding  the  nonlinear  shock  discontinuity  and  underlaid  by  a thin 
viscous  disturbance  sublayer  as  schematically  illustrated  in  Fig.  2. 

To  achieve  an  analytical  solution  in  the  leading  approximation,  the  fol- 
lowing assumptions  are  introduced  (see  Ref.  5,8  for  more  detailed  dis- 
cussion of  their  validity). 

(a)  The  incident  shock  is  a discontinuity  across  which  Rankine-Hugoniot 
shock  jump  relations  are  satisfied. 

(b)  The  nonlinear  transonic  terms  in  the  outer  inviscid  flow  regions  are 
neglected,  since  most  of  the  significant  nonlinear  transonic  effect  is  al- 
ready accounted  for  in  the  shock  jump  relations.  We  thus  deal  with  a 
linearized  (but  rotational)  boundary  value  problem  surrounding  the  non- 
linear shock  discontinuity. 

g 

(c)  Following  Lighthill’ s treatment  of  the  oblique  shock  case  , we  im- 
pose the  incident  normal  shock  jump  conditions  only  at  the  boundary 
layer  edge  and  neglect  the  details  of  the  shock  pentration  into  the 
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underlying  region;  since  the  correct  shock  pressure  jump  at  the  edge  is 
accounted  for  while  below  the  sonic  level  in  the  boundary  layer  no  dis- 
continuity can  exist,  the  shock  decay  across  the  supersonic  non-uniform 
flow  region  is  in  fact  roughly  simulated  by  this  approximation. 

(d)  The  viscous  disturbance  sublayer  is  assumed  thin  enough  to  lie 

9 

within  the  linear  portion  of  the  undisturbed  boundary  layer  profile 
Provided  is  not  so  close  to  unity  that  the  incident  shock  thickness 
becomes  a significant  fraction  of  the  boundary  layer  thickness,  this 
approximate  flow  model  contains  all  the  essential  global  features  of  the 
mixed  transonic  character  of  the  non-separating  normal  shock-turbulent 
boundary  layer  interaction  problem  including  lateral  pressure  gradient 
effects  over  a wide  range  of  Reynolds  numbers.  Moreover,  the  linea- 
rized theory  involved  is  now  amenable  to  analytical  treatment  by  ob- 
taining solutions  in  each  of  the  three  regions  shown  in  Fig.  2 and 
appropriate  matching  of  them. 

The  linearized  small  disturbance  supersonic  region  1 is  governed  by 
the  wave  equation.  Thus  if  77  = tj^x  ) is  the  small  vertical  displacement 
of  the  interface  due  to  shock  interaction  and  incoming  wave  disturbances 
are  ruled  out,  the  pressure  perturbation  about  the  local  undisturbed 
value  far  upstream  is  given  by 

(1) 

where  7^  = tj^(x  - The  subsonic  disturbance  flow  in  quadrant  3 is 

caused  by  the  interaction-generated  interface  displacement  rj^M  along 
y = 6 plus  the  post- shock  perturbations  along  x = 0+  resulting  from  the 

it 

impingement  of  region  1 Mach  wave  disturbances  on  the  shock.  The  so- 
lution for  the  local  post-shock  pressure  perturbation  p^'  - p - pq3 
which  vanishes  both  at  x = 0,  y — » 00  and  x— ♦ ao,  y = 0 can  be  found 
by  Fourier  Sine  transformation  to  be^ 


p3'  U.y) 

Po3 
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= |-  J P3(y,k)  sin  kx  dk 


(2A) 


* For  example,  the  shock  jump  relations  give  p'  3(0,y)  - p'  ^O.y) 
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where  y = y - 6 and 
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Within  the  rotational  inviscid  disturbance  boundary  layer  region  2,  the 
Fourier  transform  iT(k,y)  of  the  non-dimensional  disturbance  pressure 
p^/p^Cx.y)  is  governed  by  the  equation9 
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drj  ' drj 


k2(l  - 


M ) ij 
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where  y/6  = 77.  The  solution  obeys  the  inner  boundary  condition 
(d7r/drj)w  = 0 along  the  effective  wall  position  (see  below)  and 

the  following  outer  condition  obtained  by  pressure  matching  along  the 
interface: 

^P/p0i  Po3 

(k,  T)  = 1)  = H2(k)  = H1  (k)  + — — + ^H3  (4) 

Pol 


where  Ap  = P „ - p is  the  basic  normal  shock  jump  and  H0(k)  is  the 

Oo  O I o ^ 

Fourier  transform  of  h (x)  = p!(x,  l)/p  while  H,(k)  =J  h,e  "l  x dx, 

OO  . i £ ^ Ol  1 - 00  1 

/*•  1KX 

h3e  dx.  Then  introducing  the  "unit"  solution  Q(k<TJ)  of 
Eq.  (3)  satisfying  Q(k,0)  = 1,  dQ/drj  (k,0)  = 0,  the  general  solution  with 
the  stated  boundary  conditions  in  physical  variables  is 
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Within  the  viscous  disturbance  sublayer,  a well-known  solution3  for  a 

linear  basic  flow  profile  U =>*  (t  j\i  )y  yields  the  sublayer  displacement 

o w w 

thickness  (effective  wall  position): 
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where  x^  is  the  large  scale  upstream  influence  parameter  defined  by 
Lighthill  in  terms  of  the  properties  of  the  function  n<x^y). 


Applying  these  regional  solutions  to  the  interfacial  matching  conditions 
of  continuous  pressure  and  streamline  slope  at  the  boundary  layer  edge 
yields  integro-differential  equations  for  the  transform  of  the  interaction 
pressure  along  the  interface  both  upstream  and  downstream  of  the  shock; 
the  details  including  their  lengthy  but  straightforward  solution  by  oper- 
ational methods  are  described  in  Refs.  5,  8,  Inversion  by  the  method 
of  residues  and  use  of  the  convolution  theorem  then  yields  the  physical 
pressure  perturbation  distributions  along  the  interface  and  also  along 
the  wall.  In  particular,  the  wall  disturbance  pressure  is  found  to  have 
the  following  asymptotic  behavior  far  upstream  and  downstream,  re- 
spectively: 
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Eq.  7 A describes  an  exponential  pressure  decay  with  upstream  distance 
with  a characteristic  upstream  distance  scale  ~(*J  whereas  (7B) 
predicts  far  downstream  an  algebraic  (1/x)-  type  decay  of  the  local 
disturbance . 

Under  the  present  assumptions,  Inger  s solution  for  the  shear  stress 
perturbation  within  the  viscous  disturbance  sublayer  can  be  carried  over 
directly  to  the  present  problem;  in  terms  of  the  surface  pressure 
and  a function  F(x)  of  the  properties  of  n(K,y)  as  follows- 
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Although  this  integral  has  a branch  point  at  the  origin  whic  = makes  it 
very  difficult  to  evaluate  for  the  upper  complex  plane  values  required 
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for  a downstream  solution  ( x > 0),  the  integration  can  be  done  by 

residues  for  the  important  region  under  and  upstream  of  the  shock 

with  the  result  expressed  in  terms  of  p ’ . 

w 

2.2  Typical  Features  of  the  Zero  Mass  Transfer  Solution 

The  unit  solution  it  r Q(k,r})  of  Eq.  (3)  was  obtained  by  outward  numer- 
ical integration  using  a Mach  number  profile  MQ(y)  based  on  the 
accurate  turbulent  model  of  Inger  and  Williams^  • which  was  developed 
especially  for  this  purpose  and  is  in  good  agreement  with  transonic  ex- 
perimental data  (Fig.  3).  It  can  be  further  extended  to  include  the 
presence  of  surface  mass  transfer  as  shown  below. 

As  originally  constructed,  the  foregoing  theory  was  concerned  with 
establishing  a sound  analytical  structure  and  a reliable  numerical  im- 
plementation without  regard  to  computational  optimization.  However,  by 
subsequent  streamlining  it  has  been  possible  to  reduce  the  calculation 
time  enormously.  For  example, the  use  of  only  two  terms  in  evaluating 
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the  residue  series  summations  involved  yields  negligible  (.  1 $ or  less) 
error  in  the  pressure  for  > 1.05  while  cutting  the  time  by  two- 
thirds.  Furthermore,  use  of  an  approximate  closed  form  expression  for 
the  Mach  number  profile  that  provides  a direct  analytical  relation  for 

the  profile  slope  dM_/dy  without  numerical  differentiation  further  im- 

° 12 
proves  the  computational  efficiency  by  an  order  of  magnitude  Such 

improvements  result  in  an  IBM  370-158  running  time  of  10-15  seconds 
for  one  complete  interaction  case  including  all  pressure  distributions, 
boundary  layer  thickening  and  skin  friction,  which  provides  an  eco- 
nomical research  tool  for  parametric  studies  and  applications. 

Typical  features  of  the  resulting  interaction  pressure  field  solution  are 
illustrated  in  Fig.  4,  where  both  the  boundary  layer  edge  and  wall 
pressure  distributions  are  shown.  The  streamwise  extent  of  the  inter- 
action is  seen  to  contract  noticeably  with  increasing  shock  strength, 
accompanied  as  expected  by  a strengthening  of  the  local  pressure  jump 

at  the  boundary  layer  edge;  these  trends  agree  with  experimental  fin- 
7 1 3 

dings  .Moreover,  it  is  seen  that  the  shock-induced  lateral  pressure 
gradients  are  significant  within  a region  of  several  boundary  layer  thick- 
nesses upstream  and  downstream  of  the  shock  foot,  the  wall  pressure 
being  higher  then  the  edge  ahead  of  the  shock  and  lower  behind  it. 

Further  upstream  where  the  wall  and  edge  pressure  equalize  they  decay 
exponentially  with  distance.  Along  the  boundary  layer  edge,  a local 
pressure  jump  occurs  across  the  shock  at  x = 0 followed  by  a small 

region  of  subsonic  post-shock  expansion  and  subsequent  recompression, 

7 13 

again  in  qualitative  agreement  with  experimental  observations  ' . Far 

downstream,  the  wall  and  edge  pressures  again  equalize  and  rise  mo- 
notonically  to  the  final  post- shock  level  like  1/x. 

It  is  noted  that  the  local  shock  jump  at  the  boundary  layer  edge  and 
its  rapid  lateral  smoothing  across  the  underlying  subsonic  flow  region 
that  yields  a continuous  wall  pressure  distribution  are  important  physical 
features  that  cannot  be  accounted  for  without  considering  the  dp/dy  effect. 
Moreover,  inclusion  of  lateral  pressure  gradients  captures  another  inte- 
resting feature:  the  existence  of  a subsonic  post-shock  expansion  region 


at  the  boundary  layer  edge  (here  due  to  the  change  across  the  normal 
shock  of  the  upstream  compression  waves  from  the  interaction-induced 

boundary  layer  thickening,  see  Fig.  1).  This  is  qualitatively  confirmed 
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by  both  experiments  and  detailed  numerical  solutions. 

To  illustrate  the  accuracy  and  limitations  of  the  present  theory.  Fig.  5 

shows  a comparison  between  its  wall  pressure  predictions  and  a more 

14 

elaborate  numerical  solution  and  the  pipe  flow  experimental  data  of 

i n g 

Gadd  for  an  adiabatic  flow  at  M = 1.  12  and  Re^  = 9.6  x 10  . It  is 
seen  that  there  is  good  agreement  between  the  two  theories  when  un- 
corrected for  the  post-shock  channel  blockage  effect  due  to  the  inter- 
action-induced  boundary  layer  thickening;  each  possesses  virtually  the 
same  accuracy  compared  to  experiment,  both  moderately  underestimat- 
ing the  upstream  influence  and  overestimating  the  rate  of  post-shock 
pressure  recovery  as  is  typical  of  all  such  channel  flow  experiments 
Thus,  under  exactly  the  same  conditions  the  present  theory  predicts 
the  interactive  wall  pressure  with  essentially  the  same  accuracy  as 
does  a far  more  elaborate  and  expensive  numerical  solution,  while  also 
accounting  for  the  subsonic  post-shock  expansion  at  the  boundary  layer 
edge;  hence  it  provides  a sound  basis  for  interpreting  experimental  data 
and  further  extension  for  unseparated  flows. 

2.3  Effect  of  Suction  or  Blowing 

Some  preliminary  study  discloses  that  surface  mass  transfer  can  in- 
fluence the  interaction  in  numerous  ways,  which  can  be  conveniently 
broken  down  by  regions  as  follows. 

1. )  Basic  Undisturbed  Flow 

a.  Alter  values  of  t ,6 

wo  o 

b.  Change  damping  of  turbulence  across  laminar  sublayer 

c.  Change  profile  shape  away  from  wall 

2.  ) Viscous  Disturbance  Sublayer  Solution 

New  mass  transfer  - induced  profile  curvature  and  normal 
velocity  terms  in  disturbance  equations 

3.  ) Pressure  Field  in  Inviscid  Rotational  Disturbance  Region 
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New  non-parallel  mass-transfer  terms  in  governing 
disturbance  equation 

The  effects  under  1.)  are  by  far  the  most  important  and  will  be  dealt 

with  below.  As  regards  the  others,  they  may  be  neglected  under  the 

assumed  conditions  of  small-to-moderate  normal  mass  transfer  rates 
. - 3 

(B  = m /p  U < 10  ) typical  of  practical  applications,  according  to 

w ef  e( 

the  following  considerations.  Under  the  continued  assumption  that  the 

viscous  disturbance  sublayer  lies  within  the  laminar  sublayer  region  of 

a turbulent  boundary  layer,  the  mass  transfer  effect  in  the  leading 

approximation  does  not  alter  the  U0(y)  profile  curvature  but  only  its 

slope  (r  );  moreover,  this  slope  change  is  well-approximated  as  a 

linear  function  of  B.  Furthermore,  studies  of  solutions  to  the  hydro- 

15 

dynamic  stability  equations  for  large  Reynolds  numbers  have  shown 
that  including  the  Vq  term  for  small  to  moderate  B has  only  a very 
weak  effect  on  solutions  for  "parallel  shear  flow"-type  problems;  since 
our  equations  are  very  similar  in  form  we  infer  that  the  explicit  Vq 
terms  may  also  be  neglected  in  solving  both  the  viscous  disturbance 
sublayer  and  the  overlying  inviscid  pressure  perturbation  equations. 

Thus,  to  a consistant  order  of  approximation  the  form  of  these  equa- 
tions is  not  significantly  changed  by  mass  transfer  but  only  the  input 
values  from  the  incoming  basic  boundary  layer  flow,  provided  B is  well 
below  the  so-called  "blow-off"value. 

Consider  now  the  modeling  of  mass  transfer  effects  on  the  incoming 

flow,  assuming  for  simplicity  that  any  such  blowing  or  suction  is  on  the 

average  uniform  and  normal  to  the  wall  and  that  its  streamwise  extent 

is  large  compared  to  the  short  interaction  range,  extending  far  enough 

upstream  to  have  established  a well-defined  local  equilibrium  profile 

in  the  incoming  boundary  layer.  Then  the  mass  transfer  effect  on  skin 
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friction  can  be  reliably  described  by  the  relation1 
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is  the  zero-blowing  value*.  According  to  Eq.  (9),  suction  for  example 
increases  t and  this  should  have  an  influence  on  the  interaction 
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that  is  qualitatively  similar  to  that  of  decreasing  the  Reynolds  number. 
The  associated  mass  transfer  effect  on  boundary  layer  thickness  can  be 
calculated  as  follows.  Since  it  is  well-known  that  the  momentum  thick- 
ness to  boundary  layer  thickness  ratio  and  the  Crocco  energy  equation 
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solution  are  both  insensitive  to  moderate  amounts  of  mass  transfer,  we 
can  use  the  following  approximate  zero  blowing  relationship11  based 
on  a power-law  (IJ  ~ y 1 ) profile: 
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where  0 * is  found  from  the  momentum  equation  inclu- 

ding mass  transfer1  ‘ : 
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and  X here  is  the  running  length  from  some  upstream  reference  point. 

Then  under  the  aforementioned  assumption  that  B is  a constant  over 

some  region  X < X<  X 2 (and  zero  outside)  with  X^,  X ^ far  in  front  of 

and  behind,  respectively,  the  interaction  zone  X =“  L,  and  using  Eq.  (9) 

plus  the  approximate  power  law  formula  C.  ~ X ^ to  facilitate 

to 

analytical  integration,  Eq.  (1 1 ) yields 

!_  =,  (N+3)  [ 23(N+3)  l/Tref/Te)c  ^ 1 

L 2(N+1 ) * Cfo  B l 1 ‘ 2 (N+2)  V\  2 / fo  J ' 


* 

We  have  used  for  this  a reference  temperature  - modified  Schultz- 
Grunow  relation 


with  N =*  7 to  10.  Use  of  this  result  in  Eq.  (10)  thus  enables  estimation 
of  the  mass  transfer  effect  on  boundary  layer  thickness. 

Aside  from  the  foregoing  blowing  effect  on  the  profile  wall  slope  (r  ), 

w 

there  is  an  additional  distortion  of  its  detailed  shape  away  from  the 
wall  which  derives  from  the  following  approximate  turbui_.it  boundary 
layer  shear  stress  distribution  recommended  by  Conrad  and  Donaldson^; 
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with  tj  = y/6  and  where  to  a good  first  approximation  for  weak  to  mo- 
derate mass  transfer  U =*  (U)  in  the  last  term.  Eq.(13)  is  used  in 

connection  with  the  basic  relation  defining  turbulent  shear  stress,  which 
can  be  written  in  the  form 
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where  it  remains  to  specify  the  total  kinematic  viscosity  distribution 
v in  the  turbulent  flow.  Now,  the  available  experimental  evidence 
suggests  that  the  form  of  the  turbulent  eddy  viscosity  relation  is  signi- 
ficantly affected  only  by  relatively  large  amounts  of  surface  mass  trans- 
17,  19 

fer  ’ ; provided  we  account  for  the  mass  transfer  effect  on  the  para- 

meters involved  (e.g.  , j ) we  thus  can  continue  to  use  the  existing 
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B = 0 viscosity  formulation  for  the  weak  - to  - moderate  blowing  or 
suction  cases  of  interest  here.  Thus  we  have  the  two-layer  piece-wise 
continuous  model^: 
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where  r is  the  Klebanoff  intermittency  factor  and  rj*  is 
found  by  requiring  the  integral  of  Eq.  14  to  satisfy  the  no  slip  condi- 
tion U(0)  = 0 upon  integrating  forward  from  the  outer  initial  condition 
U (1 ) = Ug.  Thus  the  substitution  of  Eqs.  (13)  and  (15)  plus  the  Crocco 
integral  for  the  temperature  profile  T(u)  into  Eq.  14  and  subsequent 
integration  yields  an  accurate  yet  fundamentally  - based  determination 
of  the  incoming  turbulent  boundary  layer  velocity  and  Mach  num- 

ber profiles  including  compressibility,  heat  transfer  and  moderate 
amounts  of  wall  suction  or  injection;  the  results  satisfy  all  the  proper 
boundary  conditions  including  vanishing  gradients  at  the  boundary  layer 
edge,  conform  to  the  Law  of  the  Wall  near  the  surface,  are  continuous 
across  the  entire  boundary  layer  with  a velocity  defect-type  behavior 
in  the  outer  part,  and  are  in  good  agreement  with  experiment  over  a 
wide  range  of  transonic  - to  - moderately  supersonic  Mach  numbers. 

2.4  Typical  Features  of  the  Interaction  with  Mass  Transfer 

Fig.  6 illustrates  the  typical  mass  transfer  effect  on  the  Mach  number 
profile,  from  which  it  can  be  seen  that  this  effect  is  dominated  by  the 
influence  of  mass  transfer  on  the  profile  shape  away  from  the  wall; 
including  only  the  effect  on  shear  stress  (wall  slope)  badly  underesti- 
mates both  the  magnitude  and  sign  of  the  profile  changes.  This  can  be 
understood  by  a study  of  Eq.  (13):  inserting  typical  values  in  the  right 
hand  side  shows  that  a given  amount  of  blowing  (for  example)  increases 
the  local  shear  stress  and  hence  velocity  gradient  far  more  than  it 
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reduces  the  wall  slope. 

We  now  examine  the  consequences  of  this  on  the  interaction  solution 
itself.  Fig.  7 gives  a typical  result  for  a M = 1.10,  Re  - 10*’  inter- 

A 

action,  showing  how  the  various  contributions  to  the  suction/blowing 
effect  influence  the  wall  pressure  distribution  (analogous  influences 
occur  on  the  displacement  thickness  and  skin  friction).  Whereas  the 
contribution  of  the  mass  transfer  effect  on  6q  is  negligable  compared 
to  that  on  r , the  effect  on  profile  shape  is  quite  large  (in  fact  com- 
pletely opposite  to  and  overwhelming  the  r -effect)  as  indeed  would  be 

w 

expected  from  the  aforementioned  influence  on  MQ(y).  This  conclusion, 

which  was  found  to  apply  over  a wide  range  of  conditions,  is  concordant 

12 

with  the  finding  of  Panaras  and  Inger  that  the  interaction  is  quite 
sensitive  to  the  turbulent  boundary  layer  profile  form  factor. 

Referring  hereafter  to  the  complete  mass  transfer  model  with  all 
three  effects  included  we  observe  that  suction,  because  of  its  predo- 
minant effect  in  decreasing  the  Mach  number  gradient  and  hence  en- 
hacing  the  profile  "fullness"  away  from  the  wall,  reduces  the  stream- 
wise  extent  and  thickening  of  the  interaction,  making  it  appear 

more  inviscid-like  in  character  with  a steeper  adverse  pressure  gra- 
dient; thus  suction  is  qualitatively  equivalent  to  an  increase  in  Reynolds 
number.  Blowing  has  the  opposite  effect,  tending  to  spread  out  the  inter- 
action pressure  field  and  increase  the  displacement  thickness.  The  cor- 
responding skin  friction  results  (see  Ref. 20  and  below)  are  consistent 
with  these  trends:  suction  increases  the  skin  friction  level  far  upstream, 
whereas  near  the  shock  it  is  increased  because  of  the  suction  - induced 
steepening  of  the  adverse  pressure  gradient  in  this  region. 
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3.  Parametric  Studies  and  Experimental  Comparisons 


3.1  Zero  Mass  Transfer 


THEORETICAL  RESULTS 


A feature  which  is  of  great  practical  importance  is  the  scaling  effect 
(Reynolds  number)  on  the  interaction  pressure  field;  the  present  pre- 
diction for  this  is  shown  in  Fig.  8.  It  is  seen  that  there  is  a signifi- 
cant Reynolds  number  effect  even  in  the  unseparated  case:  the  extent 
of  the  interaction  upstream  and  downstream  of  the  incident  shock  de- 
creases with  increasing  Reynolds  number,  tending  toward  a solution  typ 
ical  of  the  response  to  a simple  step  pressure  rise  at  very  high  Rey- 
nolds number  in  agreement  with  the  trends  of  the  available  experimental 
7 21  22 

data  ’ ’ . Moreover,  at  the  boundary  layer  edge  the  strengths  of 

the  local  shock  jump  and  post-shock  expansion  increase  and  decrease, 
respectively,  with  increasing  Reynolds  number;  at  sufficiently  high  Re 

l—j 

the  post-shock  expansion  region  becomes  very  small  and  weak  and 
hence  difficult  to  detect  experimentally. 


Another  important  aspect  of  fundamental  and  practical  interest  is  the 

characteristic  upstream  influence  distance  (here  defined  as  the  distance 

X upstream  of  the  shock  where  the  local  interaction-induced  pressure 
up 

rise  is  only  % of  the  overall  total).  The  present  theoretical  predictions 
for  this  at  various  shock  strengths  as  a function  of  Reynolds  number 
are  shown  in  Fig.  9,  plotted  as  the  ratio  of  X to  the  basic  (non-inter- 
acted)  boundary  layer  thickness  6 (which  also  of  course  experience 


Mach  and  Reynolds  effects).  These  values  are  of  order  unity(  X 


up 


~ 6 ) 
o 


as  we  should  indeed  expect  for  the  "short-range"  type  of  interactions 
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characteristic  of  turbulent  boundary  layers  . It  is  seen  from  Fig.  9 
that  the  upstream  influence  decreases  markedly  with  both  the  shock 
strength  and  Reynolds  number,  in  agreement  with  both  experimental 
observations  and  numerical  simulations  in  transonic  flows  . At 
moderate  Reynolds  numbers,  /6q  decreases  monotonically  with  Re^ 
approximately  as  a power  law,  whereas  it  tends  to  become  independent 


of  Reynolds  number  (perhaps  even  increasing 


slightly  with  Re  ) 


at  very  high  Re^  . These  conclusions  agree  with  the  predictions  of 
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strong  interaction  theory  for  turbulent  flows  as  recently  shown  by  Inger 
Using  the  pressure  solutions  in  the  y-momentum  equation  and  then  in- 

C O 

tegrating  twice  with  respect  to  x along  y = 6 yields  detailed  expressions  ’ 


for  the  interaction-induced  interface  displacement  (boundary  layer 

displacement  thickness  change)  AY/6q.  This  interaction  growth 

of  the  boundary  layer  is  often  of  practical  interest,  expecially  in  the 

interpretation  of  channel  flow  interaction  experiments  as  discussed  in 
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detail  elsewhere  . The  influence  of  Reynolds  number  on  this 

thickning  is  shown  in  Fig.  10.  It  is  seen  to  decrease  significantly 
with  increasing  Reynolds  number,  which  is  consistent  with  the  afore- 
mentioned pressure  distribution  trends  and  again  in  qualitative  agree- 
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ment  with  experiment  ' and  numerical  simulations  . It  is  to  be 

noted  that  the  downstream  asymptotic  values  are  5 to  8 times  larger 

when  expressed  in  terms  of  displacement  thickness  5 * instead  of  6 , 

o o 

so  that  even  the  linearized  theory  prediction  of  the  interaction  effect 
on  downstream  boundary  layer  displacement  thickness  is  generally  quite 
significant  (ranging  from  25$  to  50$  under  typical  experimental  con- 
ditions) . 


We  now  turn  to  the  disturbance  skin  friction  upstream  of  the  shock. 

The  results  of  a parametric  study  of  this  important  property  are  pre- 
sented in  Fig.  11,  which  shows  the  Reynolds  number  effect  on  the  local 

ratio  C (x)/C  at  various  shock  strengths  and  illustrates  how  C typically 
i tQ  t 

decreases  toward  the  shock  owing  to  the  adverse  pressure  gradient 
disturbance  induced  by  the  shock-boundary  layer  interaction.  It  is  also 
seen  that  increasing  shock  Mach  number  enhances  this  drop  in  skin 
friction  owing  to  the  stronger  local  interaction  pressure  gradient  in- 
volved. In  fact,  when  the  interaction  is  strong  enough,  the  present 
theory  predicts  vanishing  skin  friction  below  or  ahead  of  the  shock  at 
lower  Reynolds  numbers.  Although  this  result  tends  to  overpredict  the 
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interaction  effect  on  skin  friction  owing  to  the  linearized  small 
disturbance  approximation  in  Eq.  8 (see,  e.  g.  , the  discussion  in  Ref.  9) 
the  theory  is  nevertheless  useful  for  inferring  basic  qualitative  trends. 
For  example,  it  can  be  seen  that  Reynolds  number 
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has  the  expected  large  influence  on  the  skin  friction  and  incipient  se- 
paration behavior:  the  relative  effect  of  the  interaction  at  a given  shock 
strength  decreases  significantly  with  Re  while  separation  moves  notice- 

Lj 

ably  upstream  with  increasing  shock  strength  and  decreasing  Reynolds 

number,  all  of  which  are  in  qualitative  agreement  with  experimental 
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trends  ’ and  the  results  of  Navier-Stokes  numerical  calculations 

EXPERIMENTAL  COMPARISONS 

Although  an  appreciable  body  of  data  on  transonic  shock- boundary  layer 
interactions  has  accumulated,  most  of  these  tests  involve  high  Mach 
numbers  ( M ^ 1.3  - 1.4)  with  a distinct  lambda-shock  interaction  pat- 
tern and  definite  boundary  layer  separation  and  hence  cannot  be  com- 
pared meaningfully  with  the  present  theory.  However,  there  are  a few 
unseparated  experimental  cases  with  which  direct  comparisons  are  pos- 
sible. 

Inclusion  of  the  d p/d y effect  in  the  present  theory  enables  a comparison 
of  its  predictions  for  the  local  interaction-pressure  jump  across  the 
shock  at  the  boundary  layer  edge  versus  shock  strength  against  the 
various  experimental  values  tabulated  in  Ref.  26  (Fig.  12).  Within  its 
realm  of  validity  (unseparated  flow,  M^<  1.3)  the  theory  is  in  good 
agreement  with  the  data  and  shows  the  correct  trend  of  approaching  the 
full  inviscid  Rankine-Hugoniot  value  with  increasing  M^  or  Reynolds 
number. 

A careful  study  of  the  available  NAE  wind  tunnel  tests  of  supercritical 
27 

airfoil  sections  identified  two  interaction  cases  suitable  for  compa- 
rison. The  measured  pressure  distributions  and  correspondig  theoreti- 
cal predictions  (based  on  the  local  pre-shock  Mach  number 

and  Reynolds  number  conditions  at  the  experimentally-observed  shock 
location)  are  shown  in  Fig.  13.  The  theory  is  seen  to  predict  the  up- 
stream influence  wellj  whereas  it  overestimates  the  pressure  recovery 
downstream.  This  is  typical  of  such  airfoil  tests  and  is  apparently  caused 
by  the  fact  that,  in  contrast  to  the  normal  incident  shock  theoretically 
assumed,  the  actual  shock  occurring  in  airfoil  experiments  is  usually 
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oblique  (albeit  still  with  subsonic  post-shock  flow)  owing  to  the  non-uni- 
form  nature  of  the  surrounding  inviscid  flow  (further  discussion  of 
this  will  be  found  in  a forth-coming  paper);  the  actual  overall  shock 
pressure  rise  in  transonic  flow  can  thus  be  20  - 30  $ lower  than  the 
normal  shock  value  at  the  same  incoming  flow  Mach  number.  Never- 
theless, viewed  overall  the  theory  predicts  the  major  features  of  the 
interaction  fairly  well.  The  final  zero-mass  transfer  comparison  that 
can  be  shown  involves  the  wall  pressure  distributions  at  three  Mach 

numbers  from  the  classical  interaction  experiments  of  Ackeret,  Feld- 
7 

man  and  Rott  ( see  Fig.  14).  The  general  features  of  the  interaction 
are  seen  to  be  correctly  predicted;  as  expected,  the  linearized  theory 
increasingly  underestimates  the  upstream  pressures  with  increasing 
shock  strength.  The  overestimated  downstream  pressure  recovery 
observed  in  these  examples,  which  grows  worse  with  increasing  shock 
strength,  is  traceable  to  the  channel  flow  blockage  effect  of  the  inter- 
action - induced  boundary  layer  thickening,  which  reduces  the  effective 
14 

shock  strength 

3.2  Mass  Transfer  Effects 
THEORY 

The  influence  of  the  basic  mass  transfer  parameter  B on  the  various 
physical  properties  of  a typical  interaction  case  is  shown  in  Figs.  15A- 
D,  from  which  a number  of  interesting  conclusion  can  be  drawn.  Even 
moderate  amounts  of  suction  significantly  reduce  the  overall  stream- 
wise  extent  of  the  interaction  and  steepen  the  adverse  wall  pressure 
gradient,  whereas  blowing  has  equally  the  opposite  effect.  Concordant 
with  these  trends,  suction  also  strengthens  the  local  shock  jump 
at  the  boundary  layer  edge  while  reducing  (perhaps  even  elim- 
inating at  high  enough  B)  the  degree  and  extent  of  the  post-shock  ex- 
pansion region  (Fig.  15B).  The  corresponding  thickening  effect  of  the 
interaction  shown  in  Fig.  15C  is  very  significantly  influenced  by  mass 
transfer;  for  example,  the  moderate  suction  value  B = - .0003  reduces 
A Y (oa)  nearly  five-fold.  The  influence  of  B on  the  interactive  skin 
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friction  (Fig.  15D)  is  composed  of  two  opposing  effects:  far  upstream, 
the  skin  friction-increasing  effect  of  suction  dominates  (which  tends  to 
delay  separation)  whereas  closer  to  the  shock  the  suction-induced 
steepening  of  the  local  adverse  pressure  gradient  becomes  important 
tending  to  reduce  the  local  C^.  The  former  effect  tends  to  dominate  at 
small  B where  suction  significantly  inhibits  separation;  however,  if  B 
is  sufficiently  large  it  is  possible  in  some  weaker  shock  cases  that  the 
latter  effect  takes  over  and  suction  actually  has  a slightly  adverse 
effect  on  separation  (see  Fig.  16). 

It  is  noted  that  the  influence  of  Reynolds  number  on  the  mass  transfer 
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effect  was  also  studied  ; suffice  it  here  to  state  that  it  is  similar  to  the 
zero  mass  transfer  case  for  the  values  of  B cited  in  Fig.  15. 

The  mass  transfer  effect  on  the  upstream  influence  distance  for  various 
shock  strengths  and  Reynolds  numbers  is  shown  in  Fig.  17:  suction  re- 
duces this  influence  more  effectively  for  weaker  shocks,  the  effect  being 
linear  in  B for  very  small  B.  The  corresponding  downstream  boundary 
layer  thickening  effect  vs.  B is  presented  in  Fig.  18,  showing  the 
strong  influence  of  mass  transfer  on  the  overall  interaction. 

EXPERIMENTAL,  EVIDENCE 

To  the  authors  knowledge,  no  quantitative  interaction  data  are  presently 

available  pertaining  to  transonic  shock-boundary  layer  interactions  in 

the  presence  of  wall  mass  transfer.  However,  some  qualitative  obser- 
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vations  have  been  made  in  a channel  flow  which  lend  support  to  the 
foregoing  results:  see  Fig.  19.  The  Schlieren  photos  in  the  Figure 
show  (a)  the  typical  interaction  pattern  and  (b)  the  qualitative  effect  of 
an  unspecified  amount  of  suction  on  it,  including  an  evident  thinning  out 
of  the  boundary  layer  and  reducing  of  the  separation  zone  as  implied 
by  the  present  theory.  The  corresponding  wall  pressures  measurements 
on  the  wall  opposite  the  interaction  (Fig.  19B),  while  only  of 
qualitative  value  in  analyzing  the  downstream  interaction  field,  never- 
theless clearly  show  evidence  of  the  theoretically-predicted  reduction 
of  both  upstream  influence  and  downstream  interactive-thickening*  due  to 

★ This  reduced  thickening  results  in  a higher  downstream  interaction  pressure 
owing  to  the  smaller  channel  flow  blockage  14, 
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suction.  Such  reductions  were  also  observed  to  be  significant  in  experi- 
ments on  a supersonic  compressive  interaction  flow29;  furthermore, 
delay  of  the  separation  in  the  upstream  free  interaction  region  was  ob- 
tained, again  in  qualitative  agreement  with  the  suction  effect  predicted 
by  the  present  theory. 


i. 

I 
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4.  Concluding  Remarks 


It  is  felt  that  the  present  work  establishes  the  value  of  the  theoretical 
model  involved  as  a tool  for  engineering  studies  of  non-separting  trans- 
sonic normal  shock  - turbulent  boundary  layer  interactions  including  the 
important  influence  of  wall  suction  or  blowing.  Nevertheless,  it  is  worth- 
while in  conclusion  to  reemphasize  certain  limitations  which  warrant 
improvement  in  future  theoretical  studies. 

The  primary  limitation  in  the  theory  is  of  course  the  linearized  disturb- 
ance approximation  that  restricts  it  to  only  weak  incident  shocks.  The 
primary  consequences  of  this  approximation  are  two-fold:  one,  the 
distortion  of  the  incoming  boundary  layer  profile  and  thickness  is  not 
fed  back  into  the  local  disturbance  solution  as  it  proceeds  along  the 
interaction,  resulting  in  an  underestimate  of  the  interaction  pressures 
near  the  shock  (e.g.  Fig.  14  ) that  increases  with  shock  strength;  two, 
the  nonlinear  viscous  effects  on  the  flow  in  the  frictional  sublayer  are 
neglected,  leading  to  an  overestimate  of  the  skin  friction  disturbance 
and  a premature  incipient  separation  prediction  that  becomes  very  large 

with  increasing  shock  strength.  As  regards  the  former,  a recent  in- 
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vestigation  by  Panaras  and  Inger  has  suggested  an  approximate  means 
of  treating  the  nonlinear  profile-distortion  effects  in  terms  of  the  over- 
all shape  factor  change;  further  work  is  in  progress  to  incorporate 
suction  and  blowing  effects.  Concerning  the  latter,  an  improved  non- 
linear "lower  deck"  theory  is  required,  for  example  along  the  lines  of 

an  approximate  boundary  layer-type  integral  method  such  as  used  by 
30 

Tu  and  Weinbaum  . Study  of  such  an  improvement  is  underway  by  the 
senior  author  to  provide  a more  accurate  inter- 

active skin  friction  throughout  the  complete  - Re^  range  of  non- 
separated  interactive  conditions.  With  such  an  improved  nonlinear  feature, 
the  important  question  of  incipient  separation  and  its  control  by  suction 
can  then  be  accurately,  fundamentally  addressed  without  need  for  the 
present-day  empiricisms. 
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Fig.  1 Transonic  Shock  Turbulent  Boundary  Layer  Interaction  Flow 
Patterns 
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Fig.  2 

Approximate  Flow  Model  of  Non-Separating  Interactions 
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Fig.  5 Typical  Theoretical  Experimental  Comparison  for  Pipe  Flow 
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Fig.  10  Interaction  - Induced  Thickening  of  the  Boundary  Layer  (tfl  =0) 


Fig.  11  Shock  - Strength  and  Reynolds  Number  Effects  on  the  Theo 
retical  Local  Interactive  Skin  Friction 
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Fig.  12  Comparison  of  Theory  and  Experiment  for  Local  Pressure 
Jump  Across  Shock  at  the  Boundary  Layer  Edge 
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A.  Comparison  of  Predicted  Local  Interaction  Pressures  with 
NAE  Experiments  for  a Supercritical  NACA  6 4A4I0  Airfoil: 
~ .70,  ReCoo  = 8 x I06 
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B.  Comparison  of  Predicted  and  Experimental  Pressures  for 
the  NACA  64A4I0  Airfoil:  = .751,  Rec„  = 35  x I06 
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Compar  ison  of  Theoretical  and  Exper  iment  Interaction  Pressure 
Distributions  on  Supercritic  al  Airfoil  Sections 
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Fig.  14  Comparison  of  Pressures  with  the  Channel  Flow  Measurements 
of  Ackeret,  Feldman  and  Rott^ 


Fig.  15B  Mass  Transfer  Effect  on  Interaction  Pressure  Along  Edge  of  the  Boundary  Laye 


Fig.  15D  Mass  Transfer  Effect  on  the  Interactive  Skin  Friction  Distribution 
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